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ABSTRACT 

The two closest Gamma-Ray Bursts so far detected (GRBs 980425 & 060218) 
were both under-luminous, spectrally soft, long duration bursts with smooth, 
single-peaked light curves. Only of the order of 100 GRBs have measured red- 
shifts, and there are, for example, 2704 GRBs in the BATSE catalogue alone. 
It is therefore plausible that other nearby GRBs have been observed but not 
identified as relatively nearby. Here we search for statistical correlations be- 
tween BATSE long duration GRBs and galaxy samples with recession velocities 
V ^ 11, 000km s~ (z — 0.0367, « 155 Mpc) selected from two catalogues of 
nearby galaxies. We also examine the correlations using burst sub-samples re- 
stricted to those with properties similar to the two known nearby bursts. Our 
results show correlation of the entire long GRB sample to remain consistent 
with zero out to the highest radii considered whereas a sub-sample selected 
to be low fluence, spectrally soft, with smooth single-peaked light curves (177 
bursts) demonstrates increased correlation with galaxies within w 155 Mpc. The 
measured correlation (28% ± 16% of the sample) suggests that BATSE observed 
between 2 and 9 long duration GRBs per year similar to, and from within sim- 
ilar distances to GRBs 980425 and 060218. This implies an observed local rate 
density (to BATSE limits) of 700 ± 360 Gpc^V^ 1 within 155 Mpc. 
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1 INTRODUCTION 

The two closest Gamma-R ay Bursts (GRBs) de- 
tected so far are GRB 980425 (iGalama et al.lll99^) and 
GRB 060218 |Cusumano et al.ll2006l ; iMirabal fc Halpernl 
2006), both long duration GRBs of exceptionally low 
luminosity. Indeed GRB 980425 is usually taken as the 
archetypal low-luminosity, low variability, soft spec- 
trum GRB. During nine years of operation, BATSE 
(the Burst and Transient Source Experiment on board 
the Compton Gamma-Ray Observatory) detected 2704 
GRBs. As is well known, there is good evidence 
for two observed classes of GRB distinguished most 
clearly by Tgo, the time taken t o collect 90% of 
a bu r st's total gamma-r ay fluence l|Kouveliotou et ail 
Il993l ; iNorris et al] Il984l ). Long duration GRBs (L- 
GRBs, T90 > 2s) show a correlation between spec- 
tral peak energy and overall isotr opic energy release 
- the E p /Ei SO or Amati relation l|Amati et al.l I2OO2T ) 
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- where softer bursts are less isotropically energetic 
(although GRB 980425 is an outlier to this relation- 
ship in the sense of being too spectrally hard for its 
isotropic luminosity, it is still softer than the majority of 
BATSE long bursts). L- GRBs also exhibit a spectral lag- 
luminosity relationship llNorris. Marani fc Bonnel]||2000l ; 



ISalmonson fc Galamal I2OO2I ; lNorriall2002l l where there is 
anti-correlation between overall luminosity and the time 
delay between the observation of corresponding light 
curve features in different energy bands. Furthermore, 
GRBs with long spectral lags (and hence low luminos- 
ity) tend to have smoother lig ht curves with broader fea- 
tures , and lower peak fluxes (INorris. Scargle fc Bonnel! 
2001). Isotropic-equivalent peak luminosity has also been 
shown to correlate with light curve varia bility, originally 
for a relatively small number of burs ts (jReichart et al.l 
2001; Fcnimoro & Ramirez- Ruiz 2000) but recently con- 
firmed by Guidorzi et alj (2005) for a larger sample, 
though with a larger scatter in the data. Though the 
exact form and tightness of the correlation is currently 
a matter of some debate (Re ichart fc Nvsewanderll2005l ; 
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iGuidorzi et all 120061 ; iLi fc Paczvriskl l2006h . there ap- 
pears little doubt that the luminosity of long GRBs cor- 
relates with variability. 

GRBs 980425, 060218 and the next nearest L- 
GRB with a well observed superno v a component to 
date GRB 031203 JGotz et all 120031 : IProchaska et alJ 
I2004 I Watson et alj|2004l ) are comprehensively reviewed 
bv lKaneko et all (|2007T i7 The properties of these nearby 
bursts, along with the correlations discussed above, sug- 
gest that under-luminous L-GRBs are likely to be spec- 
trally soft with smooth, single-peaked light curves. Thus 
we can use these properties to select those bursts most 
likely to be of intrinsically low luminosity and there- 
fore drawn from a relatively nearby population. Fur- 
thermore, as suggested by arguments based on anal- 
yses of detec tor sensitivities a nd luminosity function 
calculations dPian et all 120061; Igoderberg etail 120061 ; 



iLiang et al.ll2007l ; ICobb et alj2006h , this sub-class of low- 
luminosity bursts may be many times more prevalent in 
the local Universe than high-luminosity bursts. 

Approximately 75% of the 2704 BATSE GRBs were 
long, but due to large localisation error regions few of 
these bursts have identified hosts or redshifts. However, if 
some originated within similar distances to GRB 980425 
and GRB 060218 then it should be possible to estimate 
this fraction statistically via their distribution on the 
sky, using the same technique we used pre viously to in- 
vestig ate the distribution of short GRBs (|Tanvir et alJ 
2005). If indeed there is an unidentified population of 
under-luminous, soft, smooth L-GRBs, then restricting 
the BATSE sample to those with properties similar to 
GRB 980425 should enhance any correlation signal. 

We therefore consider below the correlation be- 
tween BATSE long bursts and t wo galaxy catalogues : 
the PSCz galaxy redshift survey |Saunders et alj|2000h 
and th e Third Reference catalogu e of Bright Galaxies 
(RC3) (|de Vaucouleurs et alJll99lh . The PSCz is based 
on the IRAS Point Source catalogue, and is less af- 
fected by dust extinction in the galactic plane than other 
redshift surveys, making it an appropriate comparison 
dataset for the all-sky BATSE dataset. Though less com- 
plete, the RC3 is useful as a further comparison set and 
a check against any possible catalogue bias. 



2 METHODS 

The long burst sample consists of the 1437 long GRBs 
with location errors ^ 10° and measur ed fluences in 
the c urrent BATSE catalogue (4B(R) iPaciesas et alJ 
(1999), including web supplement). Our galaxy sam- 
ples are drawn from two galaxy catalogues, the IRAS 
PSCz galaxy redshift survey l|Saunders et al.l |2000| ) 
and th e Third Reference Catalogu e of Bright Galaxies 
(RC3) (|de Vaucouleurs et al.|[l99ll ). Using the measured 
heliocentric recession velocities of galaxies in these cat- 
alogues, we made samples of galaxies within concentric 
spheres of recession velocity v < 2000, 5000, 8000 and 
11000 km s" 1 (corresponding to radii of 28(z = 0.0067), 



70(0.017), 113(0.027) and 155(0.037 | Mpc respect i vely^ . 
Following the method described in iTanvir et al.l l|2005h . 
we define the statistic $ for each galaxy sample com- 
pared with the long GRBs as shown in equation [T] where 
9ij is the separation between the ith burst and the jth 
galaxy, and is the error circle of the ith burst position 
(sta tistical comb i ned w ith systematic errors from model 
1 of lBriggs et al.l (|l999h ). For each comparison set (burst 
sample with galaxy sample), $ thus provides a measure 
of the overall correlation between the bursts and galaxies. 



Bursts Galaxie 



$ = 



1 JBi. 
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exp 
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Values of $ can be calibrated and used to measure 
the percentage of bursts correlated with each galaxy sam- 
ple by simulating distributions of known correlation with 
that galaxy sample. To achieve this we generated a num- 
ber (501) of distributions of random bursts (each with the 
same number of positions as the number of bursts under 
consideration, the same positional errors, and distributed 
on the sk y according to the known BATSE sky expo- 
sure map (Hak kila et alj 2003)). From these we computed 
$o and its associated dispersion. Similarly, we generated 
501 100% correlated test distributions for each galaxy 
sample (this time with "host" positions chosen randomly 
from the galaxy sample and pseudo-bursts displaced from 
these positions assuming Gaussian probability distribu- 
tions within error circles chosen randomly from the burst 
sample) which enabled us to calculate $ioo- Calibration 
plots of $/$o versus percentage correlation with associ- 
ated errors could then be produced. 

To reiterate, the closest known bursts (GRBs 980425, 
060218 and 031203) were under-luminous and spectrally 
soft with smooth single-peaked light curves. In order 
to select a burst sub-sample with properties similar to 
these bursts, three separate selections were performed 
on the L-GRB sample, based on observed fluence, spec- 
tral softness, and overall light curve shape. First, as ar- 
gued above, since long-lag (and therefore likely under- 
luminous) bursts increase from a negligible proportion 
of BATSE brigh t bursts to ~ 50% at trigger thresh- 
old l|Norrisll2002l ). we ordered the L-GRBs by total burst 
fluence and selected the low-fluence half. Similarly, to se- 
lect for spectrally soft bursts, we split the total sample 
in half by the ratio of observed fluence in BATSE en- 
ergy channels 1 (20 - 50 keV) and 3 (100 - 300 keV). 
We did not attempt further cuts on the data based on 
fluence or spectral hardness since this would not be jus- 
tified for several reasons. Firstly, the number of bursts 
selected in finer cuts would be very small for statistical 
analysis. Secondly the low fluence bursts (for example) 
will obviously contain bursts that are dim due to being 
further away and these cannot be distinguished from the 
nearby bursts - we aimed purely to increase the fraction 
of nearby bursts in the samples as predicted by the corre- 
lations guiding the cuts. Thirdly, given the overall lack of 



1 Throughout this paper we use Ho = 71km s — 1 Mpc -1 . Due 
to individual peculiar velocities, recession velocity is not al- 
ways an exact distance proxy. 
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Figure 1. Measured correlation (expressed as the percent- 
age of bursts in each sample correlated with galaxy distribu- 
tion) versus galactic recession velocity for concentric spheres 
of galaxy samples. Data points have been separated along the 
x axis for clarity, and error bars represent la errors. 



Figure 2. Measured correlation between BATSE long du- 
ration GRBs and galaxies in concentric shells of increasing 
galactic recession velocity from the PSCz catalogue. Details 
as Figure [T1 



correlation in the sample as a whole, any measured cor- 
relation is likely to be very weak and further cuts beyond 
the median would be statistically dubious. The third se- 
lection for bursts of a smooth, single-peaked nature was 
made by visual examination of the light curves of the 
entire L-GRB sample. We emphasise that this selection 
was based on simple pre- agreed criteria to select smooth, 
single peaked curves broadly similar to the known local 
bursts, not to attempt detailed selection with respect to 
small scale variation. To further minimise any subjectiv- 
ity involved, selection was performed independently by 
two of us and then arbitrated by a third. Finally, the 
bursts common to all three selections (177) were then 
used to form a low-fluence, spectrally-soft, single-peaked 
sub- sample. 
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Figure 3. Measured correlation between the combined low- 
fluence, spectrally soft, single-peaked L-GRBs and galaxies in 
concentric spheres from the PSCz catalogue. 



3 RESULTS 

Figure [T] shows plots of correlation (expressed as the per- 
centage of bursts in each sample correlated with galaxy 
distribution) between BATSE long GRBs and concen- 
tric spheres of galaxy samples from both the PSCz and 
RC3 catalogues. The results show a high degree of consis- 
tency, confirming that the correlation measurements are 
not dependent on a chance choice of galaxy catalogue. In 
practice, each volume sample in the RC3 catalogue con- 
tained between 1.5 and 2 times the number of galaxies 
in the equivalent PSCz sample. This inevitably leads to 
increased dispersion in the values of $ measured with the 
PSCz samples compared to the RC3 samples, though we 
still consider the completeness and homogeneity of the 
PSCz sample to make it a more appropriate comparison 
set. 

As can be seen, long burst correlation (particu- 
larly with the PSCz samples) remains formally consis- 
tent with zero out to the largest radii considered of 
v ^ 11000 km s _1 (w 155 Mpc), confirming that nearby 
long GRBs are indeed rarely observed events. However, 
the most probable level of correlation increases with dis- 
tance, but this could not be reliably investigated beyond 
the radii considered due to the flux-limited nature of the 



galaxy catalogues meaning that there are just too few 
galaxies in the catalogues at larger radii to be useful. 

Turning now to the sub-samples with properties sim- 
ilar to known local bursts, Figure [2] shows correlations 
with concentric shells (as opposed to spheres) of galax- 
ies with recession velocity radii of — 2000, 2000 — 
5000, 5000-8000 and 8000-11000 kms" 1 . From this fig- 
ure it can be seen that each of the low-fluence, spectrally- 
soft, and smooth light curve sub-samples exhibits broadly 
equivalent, marginally increased correlations in the sec- 
ond and fourth shells. Furthermore, a combined set con- 
taining only those bursts (177) common to all three in- 
dividual sub-samples exhibited increased correlation of 
16% ± 8% (= 28 ± 14 bursts) in the 2000 - 5000 km s _1 
shell, and 19% ± 11% (= 34± 19 bursts) in the outermost 
shell. 

Returning to examine concentric spheres of galaxies, 
Figure [3] shows the cumulative correlation versus radius 
of the combined sub-sample of L-GRBs, where it can 
be seen that 28% ± 16% of low-fluence, spectrally-soft 
bursts with smooth single-peaked light curves are cor- 
related with galaxies within w 155 Mpc, equivalent to a 
total of 50 ±28 bursts in the 9 years of BATSE operation. 

It is interesting to ask whether our analysis is more 
sensitive to individual burst /galaxy correlations, or to 
correlation with large scale structure on the sky. In or- 
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PSCz galaxy shell 


"^lOO 


^ioo 


recession velocity 


incl. host 


excl. host 


(km s" 1 ) 


(arbitrary units) 


(arbitrary units) (%) 


^2000 


1875 


1684 (90%) 


2000-5000 


2873 


2663 (93%) 


5000-8000 


2354 


2151 (91%) 


8000-11000 


1810 


1616 (89%) 



Table 1. Values of <J>ioo for PSCz galaxy samples measured 
including and not including the host of the pseudo-bursts in 
correlated simulations 
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der to address this question, we repeated our correla- 
tion analyses but this time in calculation of $ioo, for 
each correlated pseudo-burst we removed the galaxy from 
which its position was generated. Thus $100 then mea- 
sures the correlation between pseudo-bursts and the large 
scale structure around the generating host. The resulting 
structure-based <3>ioo values for the PSCz galaxy shells 
are shown in Table [1] 

As this table shows, ^ 89% of the correlation sig- 
nal measured by $100 is generated from galaxies other 
than the specific hosts of the pseudo-bursts. It would 
therefore seem likely that correlation between the real 
BATSE bursts and galaxy distributions is mainly due 
to large scale structure on the sky rather than correla- 
tion with individual hosts. Given the size of the BATSE 
error boxes, this is certainly not surprising, but means 
that our analysis is unlikely to be sensitive to individ- 
ual galaxy properties. Indeed, analyses performed weight- 
ing $ by estimated individual galactic mass or star 
fo rmation rate (estim ated using the methods described 
m Lcvan et al.l(|2006l )). showed no significant change from 
the unweighted results. However, given the often large 
error boxes of individual bursts, sensitivity to structure 
is actually an advantage. Our technique effectively mea- 
sures the correlation between two 2-dimensional signals, 
and is statistical in its nature - the signal is produced 
from the sums of contributions of many burst-galaxy can- 
didate pairs and thus enhances any inherently weak in- 
dividual correlation measurements at the expense of sen- 
sitivity to individual pair properties. This also precludes 
our ability to identify individual burst /galaxy associa- 
tions. It is also known that the ty pical hosts of L-GRB s 
are blue, faint, irregular galaxies (|Fruchter et al.l [20061 ) 
and thus many potential individual hosts may be miss- 
ing from the galaxy catalogues. 



4 DISCUSSION 

4.1 Two populations of long bursts? 

Our results confirm that nearby long GRBs as a whole 
are indeed rare events, with correlation remaining con- 
sistent with zero out to ~ 155 Mpc, but with a 10% 
upper limit (lcr), equivalent to 144 bursts. Restricting 
the L-GRB sample to those with properties similar to 
known nearby bursts increased the measured correlation 
to 28%±16% of the sub-sample (= 50±28 bursts) within 
the same radius. It is worth emphasising that this means 



Figure 4. The observed burst rate, p b s , within the enclosed 
volume out to redshift z cnc . The filled circle shows the ob- 
served burst rate within our investigated total volume. The 
area enclosed by the dashed lines represents the upper and 
lower limits of the predicted observed rate for Low-Luminosity 
GRBs as a f unctio n of enclosing volume (after Figure 5(b) of 
iLiane et al] 1120071) 1. 



the 177 burst sub-sample contains a lower limit (lcr) of 22 
bursts correlated with local large scale structure, almost 
one quarter more than the lcr upper limit (18) expected 
from the correlation rate of the sample as a whole. The 
local rate density of under-luminous L-GRBs implied by 
our result is in agreement with those calculated via de- 
tector se nsitivity and luminosity function arguments. For 
example ISoderberg et al. I i|2006n argue that sub-energetic 
bursts are 10 times more abundant than typical bright 
GRBs based on the sensitivities of BeppoSAX, HETE- 
2 and Swift to GRB 980 425 and GRB 060218, and sim- 
ilarly IPian et al.l |2006l ) (using BATSE, HETE-2 and 
Swift sensitivities) found a local rate density at least 
100 times greater than that estimated f rom c osmologi- 
cal bursts alone. In addition ICobb et al.l (|2006T ) estimate 
an event ratio ~ 10 2 between low and high luminos- 
ity bursts based on assuming the Swift population of 
high-luminosity bur sts to be complete to its mean red- 
shift. Furthermore, iLiang et al.l (|2007l ) suggest that in 
order to avoid over-predicting the number of interme- 
diate luminosity GRBs observed at low redshifts, the 
high-luminosity and low-luminosity bursts must be char- 
acterised by separate luminosity functions. They choose 
models of the same functional form (smoothed broken 
power- law), but with the coefficients of each separately 
constrained to pr oduce GRB rates consi s tent w ith those 
observed. Finally, iGuetta fc Delia Valid (|2007l ) estimate 
a local GRB rate of 200 - 1800Gpc" 3 yr _1 using a lu- 
minosity function consistent with the luminosities of lo- 
cal bursts. Figure [4] shows our results plotted with re- 
spect to predicted rates of Low-Luminosity GRBs from 
ILiang et al.l l|2007l ). The results presented here can be 
seen to provide independent confirmation of these lumi- 
nosity function based predictions. 



4.2 Comparison with supernova searches 

GRB 980425 was the first G RB to be observati onally 
associated wi th a supernova (Galama et al]|l998l ). and 
GRB 060218 l|Cusumano et al] 120061 ) is the most strik- 
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ing recent example. These nearby SN associated L- 
GRBs share similar properties (low-fluence, spectrally 
sof t and smooth singl e-peaked light curves) as argued 
bv iBloom et al.l (|l998l ). and it may therefore seem that 
a search for temporal and spatial correlation between 
bursts and supernovae would provide a means of identi- 
fying individual GRB hosts, particularly at low redshift. 
However the heterogeneity and inherent incompleteness 
of SN catalogues (due for example to the difficulty of SN 
detection in dusty environments, biased surveys and mag- 
nitude limitations) make it difficult to use them to iden- 
tify GRB host galaxies. Recent observations also suggest 
that not all long GRBs are necessarily accompanied by an 
observable supernova [ ? vnbo et al.ll2006l ; iGal-Yam et al.l 
l200d : iDella Valle et alj|2006l), and of course not all Ib/c 
SNe produce GRBs (e.g. ISoderberg et all (|2006h ^. Nev- 
ertheless several groups have attempted searches for spa- 
tial and temporal coincidences o f GRBs with SNe (fo r 
example IWang fc Wheelerl (Il998l); Ikippen et all (| 19981 ): 
iNorris. Bonnell fc Watanabd ( 19991 )). Results obtained 
vary depending particularly on any restrictions imposed 
on the burst or SN properties in the analyses, and gen- 
erally the statistics remain too poor to draw any firm 
conclusions from these and other studies. In addition, it 
has been suggested recently that the host galaxies of SNe 
type Ic without GRBs may be systematically different 
fro m those with GRBs, particularly in terms of metallic- 
ity ijModiaz et al.ll2007n . 

It is therefore apparent that using SNe searches to 
identify possible nearb y GRBs is fraught with difficulty. 
ISoderberg et~aH (|2006h estimate that only ~ 3% of Ib/c 
SNe give rise to detectable low-luminosity GRBs com- 
parable to GRBs 980425 and 060218. During the lifetime 
of the BATSE experimen t, the Asiago supernova cata- 
logue (Bar bon et al.| [T999n contains only 33 events clas- 
sified as lb, Ibc or Ic within the volume considered here. 
Given further that most supernova searches, particularly 
in the 1990s, concentrated on targetting relatively bright 
galaxies, untypical of GRB hosts, we should expect that 
few, if any, GRB-SN were discovered by chance in these 
surveys. Our approach, independent of any assumed cor- 
relation other than that GRBs must occur in or near a 
host galaxy, is therefore justified and enables the plac- 
ing of limits to the number of observed nearby L-GRBs 
based on location with respect to potential host galaxy 
distributions alone. 



bursts with properties similar to those of the known lo- 
cal L-GRBs significantly increased correlation with large 
scale structure on the sky to a level of 28% ± 16% 
(= 50 ± 28 bursts) within the same radius. 

The cumulative correlation rate out to « 155 Mpc 
suggests that BATSE most likely observed between 2 and 
9 long GRBs per year similar to, and from within sim- 
ilar distances to the closest known GRBs to date. This 
implies an observed local L-GRB rate density within 155 
Mpc of 700±360 Gpc~ 3 yr" 1 (to BATSE detection limits, 
and assuming the BATSE sky exposure fraction to be of 
the order of 0.5). This is in reasonable agreement with 
230lp n Gpc _3 yr~ 1 as calculated by ISoderberg et al.l 
(2006) via combined analyses of the sensitivities of Bep- 
poSAX, HETE II, Swift. It is also in very good agreement 
with observed rates of Low-Luminosit y GRBs predicted 
from luminosity function a rguments l|Liang et alj 120071 ; 
iGuetta fc Delia Valldl2007h . 

At the time of writing (March 2007) w 50 Swift- 
detected long duration GRBs have measured redshifts 
out of a total of « 180 detections. The mean redshift 
of an unbiased sample of 16 Swift long GRBs mea - 
sured in late 2005 was z — 2.8 (Jakobsson et al.ll2006l). 
though the highest so far is z = 6.295 ( Haislip et al.l 
l2006l : lKawai et alj2006j ) and the lowest GRB 060218 with 
z = 0.0331. Comparison between datasets from differ- 
ent missions and detec t ors is difficu lt (see for exam- 
ple iBerger et af] l|2005l ); iBandl i|200fj )'). However if we 
naively assume the Swift GRB catalogue to be similar 
to the BATSE catalogue, then in a sample of 50 bursts 
we would expect to find between and 5 within 155Mpc, 
as is the case. It may be argued that since we are look- 
ing for nearby bursts, these would be more likely to have 
measured redshifts than the bursts in general, though 
all three nearest GRBs were significantly underluminous 
events. Even so, comparing to the entire Swift sample 
of 180 long GRBs, our results suggest it contains 9 ± 9 
bursts with redshift less than 0.0367. So far, only one 
has been identified as occurring within our considered 
volume: GRB 060218. This is certainly consistent with 
our limits (particularly given the relatively small num- 
bers involved), but we may expect more low-redshift long 
bursts to be identified in the future, particularly among 
spectrally-soft, under-luminous bursts with little variabil- 
ity. 



5 SUMMARY AND CONCLUSIONS 

We have analysed the correlations of the entire BATSE 
catalogue of long GRBs with measured fluences localised 
to better than 10° (1437 bursts) with galaxy samples out 
to ~ 155 Mpc from two independent galaxy catalogues. 
We find that correlation between the L-GRB set as a 
whole and samples from both galaxy catalogues remains 
within 1 standard deviation (<r) of zero out to the high- 
est radii considered. However, selecting a sub-sample of 
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